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Multifragment emission following 129 Xe+ 197 Au collisions 
at 30, 40, 50 and 60 AMeV has been studied with multide- 
tector systems covering nearly 4n in solid angle. The cor- 
relations of both the intermediate mass fragment and light 
charged particle multiplicities with the transverse energy are 
explored. A comparison is made with results from a similar 
system, 136 Xe+ 209 Bi at 28 AMeV. The experimental trends 
are compared to statistical model predictions. 



I. INTRODUCTION 

Nuclear multifragmcntation is arguably the most com- 
plex nuclear reaction, involving both collective and inter- 
nal degrees of freedom to an extent unmatched even by 
fission. As in fission, multifragmentation is expected to 
present a mix of statistical and dynamical features. 

A substantial body of evidence has been presented in 
favor of the statistical nature of several features such 
as fragment multiplicities PHlO|], charge distributions 
pLOl, and angular distributions |ll|. Recently how- 
ever, evidence has been put forth for the lack of sta- 
tistical competition between intermediate mass fragment 
(IMF) emission and light charged particle (LCP) emis- 
sion. More specifically, it has been shown that for the 
reaction 136 Xe+ 209 Bi at 28 AMeV: a) LCP emission sat- 
urates with increasing number of emitted IMFs Jl4| ; b) 
with increasing transverse energy (Et), the contribution 
of the LCPs to E t saturates while that of the IMFs be- 
comes dominant [|l5|j ; c) there is a strong anti-correlation 
of the leading fragment kinetic energy with the number of 
IMFs emitted ]i~4|] . This body of evidence seems to sug- 
gest that beyond a certain amount of energy deposition 
most, if not all, of the energy goes into IMF production 
rather than into LCP emission in a manner inconsistent 
with statistical competition. 



Given the importance of these results in showing a 
potential failure of the statistical picture and a possible 
novel dynamical mechanism of IMF production, we have 
applied the same analysis to a set of systematic measure- 
ments of 129 Xe+ 197 Au at several bombarding energies. 
In what follows we report on 1) new experimental data 
that confirm the general nature of the observations in 
JTij ; 2) new experimental data which show trends that 
are different from those observed in 3) the effective- 
ness of gating on IMF multiplicity (Aimf) as an event- 
selection strategy; and 4) the reproduction of key results 
with statistical model calculations. 



II. EXPERIMENTAL SETUP 

LCP and IMF yields and their correlations with, and 
contributions to E t were determined for the reaction 
129 Xe+ 197 Au at 30, 40, 50, and 60 AMeV. The exper- 
iments were performed at the National Superconduct- 
ing Cyclotron Laboratory at Michigan State University 
(MSU). Beams of 129 Xe, at intensities of about 10 7 parti- 
cles per second, irradiated gold targets of approximately 
1 mg/cm 2 . The beam was delivered to the 92 inch scat- 
tering chamber with a typical beam spot diameter of 2-3 
mm. 

For the bombarding energies of 40, 50 and 60 AMeV, 
LCPs and IMFs emitted at laboratory angles of 16°- 
160° were detected using the MSU Miniball @. As 
configured for this experiment, the Miniball consisted 
of 171 fast plastic (40 /mi)-CsI(2 cm) phoswich detec- 
tors, with a solid angle coverage of approximately 87% 
of 4-7T. Identification thresholds for Z=3,10, and 18 frag- 
ments were ss2, 3, and 4 MeV/nucleon, respectively. Less 
energetic charged particles with energies greater than 1 
MeV/nucleon were detected in the fast plastic scintillator 
foils, but were not identified by Z value. Isotopic iden- 
tification was achieved for hydrogen and helium isotopes 
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with energies less than 75 MeV/nucleon. Energy cali- 
brations were performed using elastically scattered 12 C 
beams at forward angles and by using the punch-through 
points of the more backward detectors to normalize to ex- 
isting data . The energy calibrations are estimated to 
be accurate to about 10% at angles less than 31° and to 
about 20% for the more backward angles. 

Particles going forward (< 16°) were measured with 
the LBL forward array [Q, a high resolution Si-Si(Li)- 
plastic scintillator array. Fragments of charge Z=l-54 
were detected with high resolution using a 16-element 
Si(300 /xm)-Si(Li)(5 mm)-plastic(7.6 cm) array |L| with 
a geometrical efficiency of ~64%. Where counting statis- 
tics allowed, individual atomic numbers were resolved for 
2 = 1-54. Representative detection thresholds of Z=2, 8, 
20 and 54 fragments were approximately 6, 13, 21, and 27 
MeV/nucleon, respectively. Energy calibrations were ob- 
tained by directing 18 different beams ranging from Z=\ 
to 54 into each of the 16 detector elements. The energy 
calibration of each detector was accurate to better than 
1%, and position resolutions of ±1.5 mm were obtained. 

The complete detector system for these higher energies 
(LBL array + Miniball) subtended angles from 2°-160° 
and had a geometric acceptance «88% of 4-7T. As a pre- 
caution against secondary electrons, detectors at angles 
larger than 100° were covered with Pb-Sn foils of thick- 
ness 5.05 mg/cm 2 (this increased the detection thresholds 
for these backward detectors) . Both the Miniball and for- 
ward array were cooled and temperature stabilized. 

For the 30 AMeV data set, the forward going particles 
(6 = 8°- 23°) were measured by the MULTICS array 
p9| , a high resolution gas-Si-Si(Li)-CsI array. Detection 
thresholds were approximately 2.5 MeV/nucleon for all 
fragments (iT=l-54), and the resolution in Z was much 
better than 1 unit for Z <30. Energy calibrations were 
performed by directing 18 separate beams into each of 
the 36 telescopes. The calibration beams had energies 
of E/A=30 and 70 MeV, and ranged in mass from 12 C 
to 129 Xe. An energy resolution of better than 2% was 
obtained. Position calibrations of the Si elements of the 
MULTICS array were performed with the procedure of 
ref. (^o|. The angular resolution was estimated to be 
« 0.2°. Charged particles emitted beyond 23° were de- 
tected with the Miniball in a setup similar to the higher 
bombarding energies described above. The complete de- 
tector system covered approximately 87% of Att. 

Data were taken under two trigger conditions: at least 
two Miniball elements triggered or at least one IMF ob- 
served in the relevant forward array. 

Further details of the experimental setups can be found 
in refs. [ElLEI. 
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FIG. 1. The average LCP multiplicity (top panel), aver- 
age transverse energy of IMFs (solid symbols), and average 
Et of LCPs (open symbols, bottom panel) are plotted as a 
function of IMF multiplicity for the reaction 129 Xe+ 197 Au. 
at bombarding energies between 30 and 60 ^4MeV. 



TABLE I. Values are given for the approximate JVimf satu- 
ration value (along with the upper limit of the integrated cross 
section in percent), the average LCP multiplicity and average 
Ef CF in the saturation region of Fig. [j], and the maximum 
average IMF multiplicity for the top 5% of the Et selected 
events, for the reaction 129 Xe+ 197 Au. 
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FIG. 2. The average IMF multiplicity (solid symbols, top 
panel), average LCP multiplicity (open symbols, top panel), 
average transverse energy of IMFs (solid symbols, bottom 
panel), and average transverse energy of LCPs (open sym- 
bols, bottom panel) are plotted as a function of Et. 



III. COMPARISON WITH PREVIOUS RESULTS 

Following the procedure outlined in JT3], the average 
LCP yields were determined as a function of iViMF (which 
serves as a rough measure of impact parameter or energy 
deposition). Fig. |l] shows an example of such an analy- 
sis for the reaction 129 Xe+ 197 Au at bombarding energies 
between 30 and 60 AMeV. The average LCP multiplicity 
((■Wlcp)) does indeed saturate with increasing Nimf, as 
observed in Jl3] . However, the value to which (ATlcp) sat- 
urates ((-^LCp) max ) rises with increasing bombarding en- 
ergy and is listed in Table B. The IMF multiplicity at 
which the saturation occurs is approximately 4-5 at 30 
AMeV and rises with increasing bombarding energy to a 
value of 8-9 at 60 AMeV. 

The average LCP contribution to E t ((E^ 013 )) satu- 
rates in a bombarding energy dependent fashion as well 
(see (# t LCP ) max in Table § and open symbols of Fig. [j], 
bottom panel). In contrast, the average IMF contribu- 
tion to E t ((i?J MF )) rises linearly with increasing IMF 
multiplicity. The significance of the bombarding energy 
dependence of these observations will be discussed in the 
next section. 

We now explore the dependence of these same vari- 
ables on E t . According to the procedure outlined in [jl5| , 
the average yields of multiplicity and transverse energy 
for both IMFs and LCPs were determined as a function 
of E t (which serves as a measure of impact parameter or 
energy deposition [^3|-^5|). In Fig. ||are plotted (Nimf), 
(Nlcp), (Ej MF ), and (£ t LCP ) as a function of E t for 
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FIG. 3. The average transverse energies of IMFs (squares) 
and of LCPs (diamonds) are plotted as a function of Et for 
the reaction 136 Xe+ 209 Bi at 28 AMeV (taken from ref. p|). 



bombarding energies between 30 and 60 AMeV. All the 
observables rise with increasing Et , in disagreement with 
the observations in [ fl5| |. In the value of (£^ CP } is 
observed to saturate to a relatively small value compared 
to ^_Ej MF )(see Fig. |J), which is at variance with the ob- 
servations in Fig. |^. The origin of this disagreement will 
be discussed in the next section. 

Lastly, according to the procedure in Jl4| , the average 
kinetic energy of the projectile-like fragment ((E/A) PLF , 
defined as the heaviest forward-moving particle in an 
event, with Zplf > 10 and 9 < 23°) has been determined 
as a function of Njmf, an example of which is given in 
Fig. |]. Here, we confirm the observation in |Q. For 
increasing iVrMF, the energy per nucleoli of the leading 
fragment decreases continuously. 

The three aforementioned observations have been used 
to suggest that, above a certain excitation energy, the 
IMFs get the lion's share of the energy while the LCPs 
lose their capability to compete |lj,|l5|. In the following 
section, we explore each of these observations and suggest 
possible alternative explanations. 



IV. INTERPRETATION 

We begin with the saturation of (-/Vlcp) and (E^ CP ) as 
opposed to the continuous rise of (El MF ) observed in 
Fig. [l]. (E"( MF ) rises linearly since 



sin 



imf \e 



IMF^ 
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where (ej MF ) is the average transverse energy of an IMF. 
Thus, the reason for the continuous rise of (i?J MF ) can 
be understood quite simply. But what is the reason for 
the saturation of (-E t LCP ) and (A L cp)? We believe that 
the values of Njmf where (./Vlcp) an d (£"^ CP ) saturate 
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FIG. 4. Top panel: the average kinetic energy per nu- 
cleoli of the projectile-like fragment is plotted as a function 
of ./ViMF (solid circles) and Nlcp (open symbols). Bottom 
panel: Same as top panel but selected from events within the 
indicated range of Et- 



represent the tails of the IMF multiplicity distribution 
which are determined by the most central collisions. 

For example, the values of IMF multiplicity at which 
the observables in Fig. [j] saturate (Nf^ F ) can be under- 
stood in terms of an impact parameter scale. Consider 
the probability P of emitting N\mf and its integrated 
yield 



IMF) 



(2) 



=N I] 



as shown in Fig. |5| for the reaction y Xe+ Au at 
50 AMeV. Average impact parameter scales, as they 
are commonly employed, are proportional to ^/S [ p3[ . 
Note that the multiplicities at which saturation occurs 
represent roughly 5% of the total integrated cross sec- 
tion (dashed line in the bottom panel of Fig. ^|). The 
iViMF value Nf$ F for which S ~ 0.05 is listed in Ta- 
ble H for each of the different bombarding energies. 
A/jjJjf tracks rather well the maximum average IMF mul- 
tiplicity ((^ViMF) max ) measured for the most central col- 
lisions (top 5% of events) based upon the E t scale. 

The above observations demonstrate that large IMF 
multiplicities (AiMF>(AiMF) max ) have small probabili- 
ties and represent the extreme tails of events associated 
with the most central collisions. In other words, events 
with increasing values of N\mf in the saturation region 
of Fig. |l| do not come from increasingly more central col- 
lisions where more energy has been dissipated. Thus, 




FIG. 5. Top panel: Probability to emit N\mf from the re- 
action 129 Xe+ 197 Au at 50 AMeV. Bottom panel: Integrated 
probability to emit JVimf or more IMFs. 



A'imf is useful as a global event selector over only a very 
limited range. 

Consequently, it is expected that statistical models 
should exhibit similar trends as those observed in Fig. [I]. 
Examples of such predictions are shown in Fig. || for the 
statistical multifragmentation model SMM (open sym- 
bols) and for percolation (solid symbols) ^7|. In 
both models an excitation energy (E) distribution was 
used to mimic an impact parameter (b) weighting. As- 
suming that 6=0 events give rise to the maximum excita- 
tion energy (E max ), we have chosen the number of events 
at a given E proportional to E m&yi — E. The "excitation 
energy" for the percolation calculation is essentially rep- 
resented by the number of broken bonds and is calculated 
as per ref. p7[ . 

Both calculations show a saturation of (Alcp) when 
plotted as a function of IMF multiplicity. This behavior 
can be understood in terms of a simple model. Consider 
the statistical emission of two particle types with barriers 
B\ and B2 (and E>2 > B\). Assume the emission proba- 
bilities are pi oc exp [—Bi/T] (i = 1, 2) with p± -\-*p% = 1. 
With the temperature T characterized in terms of the to- 
tal multiplicity n to t = n\ + "2 = ctT, and ignoring mass 
conservation, the solution for (ni) as a function of 712 can 
be calculated for a distribution of excitation energies like 
that described above. The solution of this model is shown 
by the asterisk symbols in the top right panel of Fig. ^ 
for Bi=8, i?2=24, T max =10 and a = 2 (and iVi M F=«2, 
Alcp=^i). This saturation is qualitatively similar to 
that of the other statistical models listed in Fig. || and to 
the behavior observed in Fig. |]. Furthermore, the satu- 
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ration value of (Alcp), as well as the value of -/Vimf at 
which saturation occurs, both depend on the maximum 
energy used in the calculation. Consequently, for sta- 
tistical emission one expects (and observes in Fig. |l|) a 
bombarding energy dependence of the saturation which 
reflects the total energy available to the decaying system. 
These behaviors are generic features that are present in 
any statistical model p8| . 

For completeness, the IMF and LCP yields from the 
SMM calculations are plotted as a function of E t as 
well in Fig. ^| (left panels). There is no saturation 
of (£ , J jCP ) with increasing E t as was observed in [[lB). 
Instead, this calculation shows qualitatively the same 
trends as experimentally observed in Fig. 0. 

What then causes the (unconfirmed) saturation of 
(£ t LCP ) observed in 136 Xe+ 209 Bi @ (bottom panel of 
Fig. [?|)? We believe that the saturation observed in 
136 Xe+ 209 Bi is likely due to the limited dynamic range 
of the detectors used. The charged particle yields from 
the 136 Xe+ 209 Bi reaction were measured with the dwarf 
array [p9l whose thin Csl crystals (thickness of 4 mm for 
polar angle 6 = 55 - 168°, 8mm for 9 = 32- 55° and 20 
mm for 8 — 4 — 32°) are unable to stop energetic LCPs. 
For example, protons punch through 4 mm of Csl at an 
energy of 30 MeV. Consequently, their contribution to E t 
could be significantly underestimated. 

An example of the distortions that would be caused by 
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FIG. 7. The average transverse energies of IMFs (squares) 
and of LCPs (diamonds) are plotted as a function of E t 
for the reactions 129 Xe+ 197 Au at 30 AUeV (top panel), 
129 Xe+ 197 Au again but filtered with the upper energy thresh- 
olds of the dwarf array detector |^ (middle panel), and 
136 Xe+ 209 Bi at 28 AMeV (bottom panel, taken from ref. ph. 



the detector response of the dwarf array on the similar 
129 Xe+ 197 Au reaction at 30 AMeV is given in Fig. 0. 
In the top panel is plotted (£ 4 LCP ) and (£' f IMF ) as a 
function of E t as measured by the MULTICS/Miniball 
collaboration. The thicknesses of the Csl crystals from 
these detectors range from 20 to 40 mm. Protons punch 
through 20 mm of Csl with an energy of 76 MeV. In the 
middle panel of Fig. 0, the 129 Xe+ 197 Au data have been 
"filtered" using the dwarf array high energy cutoffs which 
remove high energy particles from E t . After filtering, the 
two prominent features observed in the 136 Xe+ 209 Bi data 
set jl5| (bottom panel of Fig. [?]) then appear in the fil- 
tered data. Namely, (£" t LCP ) saturates to a small value 
and (E\ MF ) becomes the "apparent" dominant carrier of 
E t . These two features are likely to be instrumental in 
origin and therefore do not warrant a physical interpre- 
tation. 

Last of all, we come to the behavior of the average ki- 
netic energy of the projectile-like fragment (E/A) PLF as 
a function of A^mf, an example of which is given in Fig. ^ 
for 129 Xe+ 197 Au at 40 AMeV (solid circles). From the 
decrease of (E/A} phF with Vimf, it was concluded that 
kinetic energy of the PLF is expended for the produc- 
tion of IMFs |14|. It was also argued that for increasing 
IMF multiplicity, the saturation of (Vlcp) represents a 
critical excitation energy value beyond which no further 
amount of relative kinetic energy between the PLF and 
TLF is converted into heat. In other words, the IMFs no 
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longer compete with the LCPs for the available energy - 
they get it all. 

One can test the consistency of this explanation by 
studying the same observable, (E/A) FLF , but now as a 
function of Ylcp (open symbols, top panel of Fig. |^). 
We observe the same dependence as that of the IMFs 
- a monotonic decrease of (E/A) p - LF with increasing 
-/Vlcp which reaches a value of f«17 MeV at the largest 
multiplicities. This behavior persists whether we re- 
strict ourselves to the saturation region (A^mf> 6, tri- 
angles) or not (open circles). The similar behavior of 
(E/A) PLF with respect to ATimf and -/Vlcp indicates that 
the LCPs do compete with the IMFs for the available en- 
ergy. 

This can be seen more clearly by pre-selecting events 
with a better global observable, E t 24 , 25|j3"o|| , as done 
in the bottom panel of Fig. |[ Once a window of E t 
is selected, a corresponding value of (E/A) pLF is also 
determined, and there is no longer any strong dependence 
of (E/A) PLF on Nimf or Ylcp- In fact, the resulting 
iViMF and -/Vlcp selections both give the same value of 
(E/A} p - LF , consistent with a scenario where both species 
compete for the same available energy. 



V. CONCLUSIONS 

In summary, we have made a systematic study of LCP 
and IMF observables as a function of IMF multiplicity 
and transverse energy for the reaction 129 Xe+ 197 Au at 
bombarding energies between 30 and 60 AMeV. 

We observe that (-/Vlcp) and (E} jCF ^ saturate as a 
function of Nimf in a bombarding energy dependent way. 
These saturations are predicted by statistical models and 
are fundamental features of statistical decay 28 1. A 
bombarding energy dependence of (./Vlcp), (E} jCF '), and 
JYj^p is expected (and experimentally observed) within 
the framework of statistical decay. 

In addition, it has been demonstrated in a model inde- 
pendent fashion that the LCPs compete with the IMFs 
for the available energy. By using E t , a more sensitive 
event selection is obtained. The analysis also demon- 
strates the limited usefulness of event classification using 
only -/Vimf- 

We do not observe a saturation of (£' t LCP )as a func- 
tion of E t at any bombarding energy. The saturation 
of (-E^ 015 ) as a function of E t observed in ref. fis] ] 
is likely due to instrumental distortions. We can ac- 
count for this saturation by filtering the present mea- 
surements of 129 Xe+ 197 Au with the experimental thresh- 
olds present in refs. |i"^ , |l5| ]. The resulting distortions 
to the data are large and induce qualitative changes in 
the trends of the data, causing an unphysical saturation 
of (-E^ CP ). Therefore, the observations listed in [ p~4|JT5| 
do not demonstrate any measurable failure of statistical 



models that would justify invoking dynamical IMF pro- 
duction by default. While the IMFs may indeed be pro- 
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duced dynamically, the observations listed in refs. 
do not provide evidence for such a conclusion. 
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